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C$ ' Abstract 

^ \ Cosmological evolution of equilibrium plasma with a condensate of U(l)-charged 

■ bosonic field is considered. It is shown that the evaporation of the condensate is 
CO . very much different from naive expectations, discussed in the literature, as well as 

from evaporation of non-equilibrium neutral condensate. The charged condensate 

■ evaporates much slower than the decay of the corresponding bosons. The evapora- 
. tion rate is close to that of the cosmological expansion. The plasma temperature, in 

contrast, drops much faster than usually, namely as the third power of the cosmo- 
iy~j ' logical scale factor. As a result the universe becomes very cold and the cosmological 

charge asymmetry reaches a huge value. 

o ; 

^! 1 Introduction 

Bosonic condensates probably existed in the early universe and played an important 
,-C< ■ role in the cosmological history. Well known examples are the classical real inflaton field, 

>• . $ [1 , or complex field, \i describing supersymmetric bosonic condensate, which carries 

baryonic, leptonic, or some other £/(l)-charge [2|. Evaporation of the inflaton produced 
particles creating primeval plasma, while evaporation of x could generate baryon or lepton 
asymmetry of the universe. Though evaporations of a real field condensate and a complex 
one share some similarities, the rates of the processes are very much different. In the case 
of the inflaton the rate of evaporation is determined by the particle production rate and 
may be quite large, while the evaporation of a charged condensate with a large charge 
asymmetry is much slower and is determined by the universe expansion rate H = a/a 
during most of its history. The impact of a large charge asymmetry on the process of 
evaporation has been considered in ref. j3] (see also [1]). A proper account of thermal 
equilibrium with a large charge asymmetry strongly changes results of refs. where 
evaporation of bosonic condensate was considered. These works are applicable to the 
case of evaporation of uncharged condensate (e.g. inflaton) but evaporation of a charged 
condensate proceeds in a much different way, determined by thermal equilibrium with a 
large chemical potential and not just by the rate of particle production into plasma. 
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Here we elaborate the approach of ref. [3] and present more accurate and detailed 
calculations. It is found in particular that the plasma temperature in the presence of 
condensate drops very fast, as T ~ 1/a 3 , where a(t) is the cosmological scale factor. After 
the condensate disappears the cooling returns to the standard relativistic law, T ~ 1/a. 
We calculated the magnitude of the cosmological charge asymmetry produced as a result 
of the evaporation of the condensate and found that it might be much larger than unity 
because of the fast cooling of the plasma. 



2 Pre-equilibrium evolution 

We consider a charged classical scalar field x a ^ the final stage of its evolution, when 
its amplitude is small enough and the potential is dominated by the quadratic term 

U{x) = m 2 \ X \\ (1) 

while quartic or higher terms can be neglected. The equation of motion for homogeneous 
X has the form: 

X + 3H X + m 2 X = 0, (2) 

where H = a/a is the Hubble parameter and a(t) is the cosmological scale factor. 

When the mass m of x becomes larger than H the solution to equation (j2J) becomes 
oscillating: 

x(t) = Xi cos m t + X2 sin mt (3) 

where the coefficients Xi,2 are slowly varying functions of time. In particular, at the 
matter dominated (MD) regime they behave as xi,2 ~ a~ 3 / 2 . Such a regime would be 
realized if cosmic energy density is dominated by the field x itself. The energy density of 
the latter is equal to 

Px = ™ 2 (IXi| 2 + IXi| 2 ) (4) 
The charge density of x is expressed through the functions xi,2 as: 

Qx = ™{\Xi?-\Xi\ 2 ) (5) 

It is clear that xi and X2 correspond respectively to the contribution of particles and 
antiparticles. For what follows we will introduce a dimensionless parameter k which 
describes the relative magnitude of the charge density: 
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(6) 



Evidently, at this stage \k\ < 1. 

The oscillating field x starts to produce lighter particles to which it is coupled (in many 
cases this can be understood as a decay of massive x-bosons at rest into light particles). 
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This process is characterized by the decay width T. Here we will not consider the kinetics 
of the decay and do not need to specify the magnitude of T. It is assumed usually, but 
often incorrectly, that the condensate decays with the rate T and to the moment when 
H < r it practically disappears transforming its energy and charge into those of light 
particles. The evolution of the energy densities of the condensate and relativistic plasma 
in this case are described by the equations: 



and indeed p c seems to disappear when the time is large enough, t ~ 1/H > 1/T. 
However, in the case when the charge density of the condensate is large, k ~ 1, a large 
chemical potential prevents from fast evaporation of a charged condensate jS| • To describe 
this phenomenon properly one needs to include the contribution of inverse reactions into 
eq. (|7J). The problem is rather complicated and we postpone it for the future. However, 
one can still use eqs. (|IJ) at the initial stage of the process when the amount of relativistic 
particles is small and inverse reactions can be neglected. 

Here we consider simpler but practically interesting case when a partial decay of 
the condensate creates thermally equilibrium plasma which cools down by the universe 
expansion. We will show that the evolution of such plasma and, in particular, the decay 
rate of the condensate is determined by the cosmological expansion rate, H, which is 
much slower than the decay rate T. The rate of the decay of a charged scalar field into 
vacuum was considered in ref. [HI]- Parameter k introduced above (JBJ) is convenient for 
the description of such equilibrium system. However, the value of k at the equilibrium 
stage should be larger than its initial value. Indeed the density of a conserved charge 
drops down as 1/a 3 , while the energy density of mixed non- relativistic and relativistic 
matter drops somewhat faster. As a result k rises and can even become larger than 1. To 
estimate its variation we use eqs. ((Zj). They can be analytically solved giving: 



where r\ = Y{t — t in ) and rj in = T in t. If the universe expansion is dominated by the non- 
relativistic condensate, the scale factor behaves as a(t) = (t/ti„) 2 ^ 3 and the integral can 
be easily evaluated. For small rj: n/ni n ~ 1 + OAi], while asymptotically n/nm ~ a, the 
latter is evident. However, as we have already mentioned above, the regime (jHJ) remains 
true only till the inverse reactions can be neglected i.e. till equilibrium is established. 

Now we will show that for a large magnitude of the condensed field, when the particle 
number density in the condensate, n c = C(t)m 3 is large, i.e. C>1, thermal equilibrium 
is quickly established when 77 1 and thus at that moment \k\ cannot noticeably exceed 
unity. We assume that initially the cosmological matter was dominated by the condensate 
of the field x an d the contribution from light particles (quarks, leptons and photons) can 
be neglected. The kinetic equation governing the production of light particles reads 



p c = -3Hp c - Tp c 

p r = -AHp r + Tp c 



(7) 





n + 3Hn Q 



TCm 3 



(9) 
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where n is the number density of the produced massless particles. The number density of 
X-particles in the condensate initially drops as C(t) = C in /a 3 and the solution to equation 
(JHJ) is easily found: 

n (t) = Tm 3 C m t -^ (10) 
a a (t) 

Thermal equilibrium would be established when rio(t) becomes of the order of its equilib- 
rium magnitude, n (t) > Uq . To this end the following condition is necessary: 

(11) 

which is evidently satisfied for a large region of the parameter values. The mass of x~ 
field is, of course, much smaller than the Planck mass. It normally can be in the range 
10 3 — 10 10 GeV. The time during which the equilibrium is established is quite short, 




r(^-to) = ^«i. (12) 

On more comment may be in order here. From kinetic equations it follows that the 
number density of relativistic particles and X"P ar ticles should decrease as 1/a 3 , if the 
system is very close to equilibrium. However, as we see in what follows, the number 
densities of quarks or leptons remains constant during rather large period, because they 
are proportional to cube of their chemical potential, eq. (|24|). which remains constant 
for large asymmetry, fi = m, where \l and m are respectively chemical potential and the 
mass of x ( see the following section). Similar behavior is true for their energy densities, 
which naively should drop down but remain constant, see eq. f|22H . This kind of evolution 
is created by deviations from equilibrium which are proportional to H/Y and small, but 
they give corrections of the order of unity to the collision integrals because of a large 
magnitude of the latter. 



(eq) 



3 Thermal equilibrium 

We consider the case when x decays into the following channel: 

X^Sq + l (13) 

where q and I are quarks and leptons respectively. Analogous process exists for antipar- 
ticles. Leptonic, Ql, and baryonic, Qb, charges are assumed to be conserved in this 
decay. Thus the accumulated charge of the condensate is transformed into charges of 
light particles. This is the essence of baryogenesis scenario of ref . [2] . 

Particle distribution functions in thermal equilibrium are described by two parameters, 
their common temperature T and chemical potentials fij, which could change with time. 
As is well known, the equilibrium distribution of fermions is: 

exp ( ^J—H±" 
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and that of bosons is 



A = exp (1ES) _ i < 15 > 

where are energies of the corresponding particles. Fermionic chemical potentials may 
have arbitrary values, while chemical potential of bosons is restricted from above by 
the value of the boson mass, fib < Tn, otherwise the distribution function would not 
be positive. If charge asymmetry of bosons is so big that even with maximum allowed 
[Mb — m such asymmetry could not be realized, then bosonic condensate must be formed 
and the equilibrium distribution becomes: 

h = (2^CS(p) + * _ (16) 

exp ^ T ) i 

where p is the three-momentum of the bosons. Thus in the case of a large charge asym- 
metry the equilibrium plasma is still described by two parameters/functions, T(t) and 
C(t), the latter being the number density of particles in the condensate in units m 3 . 

The magnitude of chemical potentials of massless fermions is determined by the con- 
ditions of thermal equilibrium and in the particular case under consideration they are: 

H q = m = fi x /A (17) 

Here we have used the equilibrium condition 

= y^/m, (18) 

in fin 

where the summation is made over all particles in initial ("in") and final ("fin") states 
participating in all relevant reactions. For the system with two conserved charges, Qb and 
Ql in the case under consideration, there should be two independent chemical potentials. 
We excluded one of them imposing the condition Q B — Q L = 0. 

The energy and, say, lepton charge density of particle species j are given respectively 

by: 

Pi = ^Jdpp'Elf^ + f^E)] (19) 

vi = ^ y d PP 2 ^ [fm+im] (20) 

Lj = f^/ dn?E[fi{iZ)-fm] (21) 

where / is the distribution function for antiparticles and Lj is the leptonic charge of 
particle j. In equilibrium chemical potentials of antiparticles are opposite to those of 
particles, p, = — /i. We assume that = Ql x = 1- If chemical potential of x is positive 
and equal to m, then (7^0 for Xi while C — for anti-%, and vice versa if fi x = —m. 
Integration over angles in eqs. (jl9H21j) is trivially performed because the distribution 
functions in homogeneous and isotropic universe are also isotropic. 
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For massless quarks and leptons the integrals can be taken analytically. The total 
energy density of all relativistic fermions plus photons is equal to: 

Prei = Pf + fH = NfT 4 + f + JL) + g T 4 = m 4 ( 7l + 72 r 2 + 7 3 r 4 ) , (22) 

where = fif/T = fi x /(AT) is dimensionless chemical potential of fermions, £ = \i x jT is 
the chemical potential of x, and 

7i = 7TTT^> 72 = 73 = 77 + 1 • ( 2 3) 



2 n 7r 2 ' ,z 2 6 ' '° 15 

Here iVy is the total number of massless fermion species, quarks and leptons included. 
For three quark-lepton families, three quark colors, and two spin states: Nf = 48 (the 
contribution of antiparticles is already included inside the brackets of eq. (J22J))- 
The charge density of massless leptons is: 

r N t T% 
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(i + ^)=- 3 (/? + Ar 2 ) 



(24) 



where = fi x /4T, Ni = 12 is the number of lepton species which include 3 families of 
particles with two spin states, and 

* - <rh- ft - s- (25) 

The subtle points are the number of right-handed neutrino states and the contribution 
of heavy t-quarks which may be absent at relatively small T created after ^-evaporation. 
In particular, if the number of neutrino species is twice smaller due to an absence of vr, 
then Ni would be 9, and Nf = 45; if instead we consider 6 neutrinos but exclude the top 
quark we would find Ni = 12 and Nf = 42. In both these cases it can be easily seen 
that the equilibrium of the reaction x — > 3g + /, together with the condition Qb = Ql, 
no longer leads us to the simple relation fJ, q — fMi — yU x /4. Indeed we find different values 
of chemical potentials of quarks and leptons in terms of that of the x particles. We have 
done calculations for these two different scenarios, with the calculated values for fi q and 
Hi, and we have found that the picture is unchanged, and all the results presented below 
are practically the same for different particle content. 

The integrals which determine the energy and charge density of massive x-field cannot 
be taken analytically but in the limit of low temperature, T <C m, their approximate 
expressions are easily found. In what follows we express everything in units of X" m ass, 
m. In particular, when condensate is non-vanishing, i.e. // = m and C ^ we obtain: 

m - % „ C + -£%L / *4 f 1 + \ l) = C + + I (26) 



m nis^-y^T=^' (27) 
m ^ * c + ^ / f^i f 1 + 1 f ) - c + a ^' 2 + ^ (28 > 

where a\ = 0.164, = 0.16, and a.3 = 8a2/T5. Notice that the condensate creates zero 
pressure because of 8(p) in the distribution function (fTEjh 
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4 Condensate decay and plasma temperature 



Since by assumption the total leptonic and baryonic charges were conserved when \ 
relaxed down to zero in ?7(l)-symmetric quadratic potential (0) and decayed into channel 
(jlHjl . each of them satisfied the conservation equation: 

Q tot + 3HQ tot = 0. (29) 

Thus leptonic or baryonic charge density evolved as: 

Qtot ~ Qin/a 3 (30) 

We have two unknown functions T(a) and C(a) if condensate is present or T(a) and 
fi(a) after condensate evaporation. To determine the law of their evolution we need 
another equation which is the law of covariant conservation of the total energy-momentum 
tensor. In Friedman- Robert son- Walker (FRW) metric this conservation law has the form: 

ptot = -3-ff (ptot + Vtot) , (31) 

where ptot includes the energy density of relativistic particles (photons and fermions), and 
the energy density of Xi either relativistic or not, depending upon the initial conditions. 
The pressure of relativistic matter is p re i = p re i/3 and the pressure of x is given by eq. 
([20)1 which in non- relativistic case turns into (|27|) . 

It is convenient to look for the evolution of C and T as functions of the scale factor 
a, taking the latter as the independent variable. In this case the Hubble parameter 
disappears from the equation and we obtain: 

Ptot = (Ptot + Ptot) /a, (32) 

where prime means derivative with respect to the scale factor a. 

The initial values of £ = m/T and C can be expressed through the initial values of 
the total energy density, p^ and, say, leptonic charge density, Q\ = Kp^ n \ As we have 
already mentioned \n\ < 1 or it may be slightly larger than 1. The values of and 

T iin) 

are presented in figures d to El as functions of p^ and k. In particular, figs. Q to 
H2 show how T in and Ci n depend upon k when p in is fixed, while figs. |U and El deal with 
their dependence upon the total pi n , with fixed k. If k ~ 1 the initial temperature is low 
and ^> 1. In particular, for k = 1 the initial temperature does not depend upon the 
initial value of p and is equal to: 

T in w 0.05m (33) 

However, even if k is slightly smaller than unity, T^ n ) becomes much higher and, for large 
Pi n , even relativistic initial state could be realized, i.e. m/T < 1. E.g. for k = 0.9 and 
Pin = 10 3 the initial temperature is T in ~ m. Such a large rise of temperature is explained 
by much larger energy release into plasma by XX annihilation, which is absent in the case 
of | | = 1 when antiparticles of x are initially absent. 

We have numerically integrated this equation with constraint (j30jl . The solution is 
divided into two parts: initially with non- zero and running C(a) (condensate), p = m, 
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and running T(a). Then at some stage C reaches zero and the condensate disappears. 
After that the solution is searched for C = but running //(a) and T(a). The details are 
presented in the Appendix. The calculations are very much simplified in the limit of low 
temperatures when T<m. The functions ptot and ptot in this limit are simple algebraic 
functions and not integrals over energy and the evolution equation turns into: 

£ 2a 2 r 5/2 + &r 2 + 47 3 r 4 i (u) 

£ 4a 2 £" 5 / 2 + 3/?i^- 2 + 4 73 e- 4 a 1 j 

Depending upon which term dominates in the £ dependent coefficients in the l.h.s. 
of this equation we obtain different regimes of the cosmological cooling. At relatively 
small £, when £ -4 dominates, the temperature drops as T ~ l/a. It is the usual law of 
relativistic cooling. If the first term, £ _5,/2 , dominates, T ~ l/a 2 . This is also the known 
regime for dominating non-relativistic matter. However, if the term £ -2 is dominant, the 
temperature drops unusually fast: 

T ~ l/a 3 (35) 

We have found that this regime is the typical one in the course of the condensate evap- 
oration, so the temperature becomes very low at the end. It is quite unusual situation. 
Normally one would expect that a source of energy created by the evaporating condensate 
would lead to plasma heating and a slower decrease of the temperature in the course of 
expansion. We observe an opposite picture. Probably this phenomenon is induced by a 
very slow decay of the condensate and efficient cooling of the particles in the plasma by 
the scattering on the "refrigerator" created by a large number of cold x-particles at rest. 

Here we would like to remark that the derivation of the temperature's evolution law 
can be extended to the case of a different particle content, even if eq. (|17|) is no longer 
true, as we have already pointed out. In fact in such a system, when both energy and 
charge densities are dominated by terms, we can generalize eq. (jHSj) to the following: 

T ~ l/a 32 (36) 

where: 

= 2N l tf + 2N f ^ f -N lfM 
Z 3N l ^ + 3N ff i 2 f -2N l p l [ ' 

If we now insert the chemical potentials given by the equations of equilibrium with respect 
to reaction ([TBI , together with the equality of charges Qb — Ql = 0, we will find for a 
different particle contents, that z — 1, which brings back eq. (|35|). 

In non-relativistic well as in the relativistic one, with dominating condensate, 

the behavior of C(a) can be predicted analytically from evolution equations (ffijj) and (|31j): 

C' = ~(C-u(T)) (38) 

where u(T) contains the temperature-dependent terms which contribute to the energy or 
to the charge density. Thus if the condensate dominates the energy and charge densities, 
i.e. C ^> u(T), the evaporation is described by the law C ~ l/a 3 . 
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The evolution of the interesting quantities is presented in figures El and [7| In accordance 
with the arguments presented above the amplitude of the condensate follows the simple 
law: 

C ~ 1/a 3 (39) 

which holds through most of its life. Only when the energy density of the condensate 
becomes approximately equal to that of the plasma, C(t) starts to drop faster and at 
some moment it reaches zero. All the plots with p = m = const have been stopped at 
the precise moment of the complete evaporation. 

If the initial temperature is high, we do not have an accurate analytical approximation 
for the energy, pressure, and charge density of \ with the large chemical potential p = m. 
In this case we need to solve differential equation (}32|) with the integrals for p, p, and Q 
taken numerically. The solutions are presented in fig. |H1 Initially the temperature drops 
as T ~ 1/a and later, when it becomes smaller than m before the complete evaporation 
of the condensate, the regime changes to that discussed above in the non-relativistic case. 

5 Results and discussion 

The obtained results are presented in figs. El El El andE3 All the quantities which appears 
in these figures are normalized to a proper power of the mass of Xi namely: 

T — > T/m, p — > p/m, 
Q^Q/m 3 , p^p/m 4 

and they are expressed as functions of the scale factor ratio a/aj n . Most of these graphs 
are in the log-log scale, except for the 'Power of T' and 'Zoom of C plots, which are in 
the log-linear scale. 

In fig. El a special case k — 1 is presented. We see in the upper left box how the 
temperature evolves, while the box below (middle left), named 'Power of T', shows that 
the cooling follows the law T ~ 1/a 3 . The same law holds for the condensate, as it 
is shown in the right upper box, that is, C ~ 1/a 3 , while the middle right box called 
'Zoom of C, is a zoom of its final evolution, and shows the instant of the complete 
evaporation. The calculated cosmological charge asymmetry, (3 = Q/n^ = Q/0.24T 3 , 
is presented in the left lower box: there are lines for the condensate (solid line), decay 
products (dashed line), and relativistic light charged particles (dot-dashed line). It can be 
seen that the produced asymmetry strongly rises during the non-relativistic regime. This 
is understandable because the temperature drops very fast, T 3 ~ a~ 9 while the charge 
density drops only as Q ~ a -3 . We have found that the ratio of (3{a — 1) to (3 at C — is, 
in this case, of the order of 10 8 . The last box (lower right) shows how entropy in comoving 
volume, S = sa 3 = a 3 (p + p)/T, evolves: we see, as it was for /3, that the fast cooling 
results in S rising as a 3 because the main contribution to the energy density comes from 
the condensate, which drops only as 1/a 3 . 

In fig. [7] the same graphs are presented, but for different initial values of temperature 
and condensate which correspond to k — 0.9, and respectively to T in w 0.65m. We see 
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that all the quantities evolve similarly to the previous case presented in figure but here 
one can better see how the behavior of temperature switches from T~ 1/a to T ~ 1/a 3 
with growing scale factor. The moment when it happens is precisely the moment when 
the £ -2 -terms in the energy and charge densities start to dominate over the standard 
ones, see eq. (|H4|) . In this situation we have found the ratio of the final to initial charge 
density to be about 10 3 . 

Next, in fig. |H] the case of initial relativistic state is presented. In order to have 
relativistic initial state we had to choose k — 0.1 and a huge initial total energy density, 
Pi n ~ 10 6 m 4 . An important thing to be noticed here is that the condensate can evaporate 
completely only when the non-relativistic limit is reached, i.e. T < m. 

Finally, in fig. E3 we present evolution of temperature, chemical potential, energy 
densities, baryonic charge, and entropy, when the condensate has completely evaporated. 
In the first box, clockwise from upper left corner, we see that the temperature follows the 
standard evolution law for a relativistic plasma, T ~ 1/a. In the second box we have the 
same graph for the chemical potential, which also behaves as fi ~ 1/a. The third box 
contains: entropy in comoving volume (thick solid line); charge density, or to be more 
precise (3, contributed by relativistic (thin solid line) and non-relativistic (thin dashed 
line 1 ) species. Finally, the fourth box (lower left corner) shows the evolution and the 
ratios of the energy densities of charged relativistic particles (solid line), photons (dot- 
dashed line), and non-relativistic ^-particles (dashed line). We see that both entropy in 
comoving volume and (3 are conserved. One may wonder if a non-zero chemical potential 
could break the entropy conservation law, but it can be easily checked that this is not so 
because every term in the expression of charge density and entropy in comoving volume 
behave as 1/a 3 . This, together with the fact that the energy density of non-relativistic 
particles (as well as their charge density) is orders of magnitude smaller than that of the 
relativistic ones, as one can see in the 4th box, gives the result just obtained. Hence the 
Affleck and Dine scenario of baryogenesis could indeed lead to a huge baryon asymmetry 
of the universe. 

To dilute the asymmetry to the observed level a subsequent release of entropy is 
necessary or evaporation into initially hot plasma. Both such mechanisms were discussed 
recently in ref. where one can find references to many earlier papers. 



In this appendix we present some details of the numerical solution of the differential 
equations 2 : 



1 This line, as well as the thin dashed one in the fourth box, seems to be vertical but in fact it is not: 
it is an extremely fast decrease of the contribution of non-relativistic x an d X particles. 

2 All the quantities here are normalized to proper powers of x-mass, e.g. T — » T/m, p — ► p/m . 



A Appendix 




(40) 
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We have to insert into these expressions the explicit values of Q tot and p tot , which depend 
either upon T(a) and C(a) when p — m, or upon T(a) and p(a) when C = 0, that is: 

Qtot = Qtot(T,C) ^ | Qtot = Qtot (T, p) 
Ptot = Ptot(T, C) \ p^t = PtotiT, p) 



In the first case, when C ^ 0, these equations read: 

(42) 



Q' tot = c + cr 
P ' tot = c + v r 

where ((T) and T)(T) are the derivatives of Q to t(T,C) and p to t(T,C) with respect to T. 
The system can be solved together with equations (|4T)|) to have: 



(77 - C) T' = 3 [(5 to i - (p tot + p tot )] /a 

(77 - C) C = 3 [C(Ptot + Ptot) - V Qtot] /a 

As for the second case, when C = 0, the derivatives of Q to t an d p tot are: 

Q' tot = f p' + gT 
P ' tot = hp' + lT' 



(43) 



(44) 



where f(p,T) and h(p,T) are the derivatives of Q tot (T,p) and p tot (T,p) with respect to 
p, while g(p,T) and f(p,T) are the derivatives of Q to t(T,p) and p tat (T, p) with respect 
to T. Again, the system can be solved together with equations (J4T)|) to have: 

(5-/1 - //) T' = 3 [/(p io t + ptot) - h Qtot] /a 

(45) 

f(gh - fl) p! = 3 [// Q to4 - gf {ptot + Ptot)] /a 

Now the system of differential equations (JH3), as well as system (jIHJ), can be solved 
numerically. 
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Figure 2: Initial values of temperature and condensate as functions of k, with the fixed value of the 
total energy density m~ 4 p to t — 100. 
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Figure 4: Initial values of temperature and condensate as functions of the energy density, with fixed 
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Figure 5: Initial values of temperature and condensate as functions of the energy density with fixed 
k = 0.1. 
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Figure 6: This figure refers to the special case k = 1. Clockwise from upper left corner: evolution 
of temperature, evolution of the condensate (upper right), the complete evaporation of the condensate 
(middle right), evolution of entropy in comoving volume (lower right), evolution of baryonic charges 
(lower left corner), and the power of T(a) (middle left), as functions of the scale factor a. The 'Baryonic 
Charges' graph has the following lines: solid line is for the condensate; dashed line is for x an d X particles; 
dot-dashed line is for relativistic charged particles. 
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Fi gure 7: This figure refers to the case k = 0.9, while m 4 ptot — 100. Clockwise from upper left corner: 
evolution of temperature, evolution of the condensate (upper right), the complete evaporation of the 
condensate (middle right), evolution of entropy in comoving volume (lower right), evolution of baryonic 
charges (lower left corner), and the power of T(a) (middle left), as functions of the scale factor a. The 
'Baryonic Charges' graph has the following lines: solid line is for the condensate; dashed line is for \ an d 
X particles; dot-dashed line is for relativistic charged particles. 
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Figure 8: This figure refers to the relativistic initial state, that is T > m. Clockwise from upper left 
corner: evolution of temperature, evolution of the condensate, the complete evaporation of the condensate, 
the power of T(a), as functions of the scale factor a. We see that the evaporation happens only when 
temperature has dropped below m. 
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Figure 9: This figure refers to the second part of the problem, when the condensate has completely 
evaporated and fi is no longer a constant. Clockwise from upper left corner: evolution of temperature, 
chemical potential, entropy in comoving volume and baryonic charges, and energy densities, as functions 
of the scale factor a. The 'Energy densities' graph has the following lines: solid line is for relativistic 
charged particles; dashed line is for photons; dot-dashed line is for x an d X particles. The 'Entropy & 
Charge' graph has the following lines: the solid line is the charge of relativistic particles; the dashed line 
is the charge of \ an d X particles, and the thick solid line is the entropy. 
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